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We study a simple sandpile model of active-absorbing state transitions in which a particle can hop 
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out of a site only if the number of particles at that site is above a certain threshold. We show that 
the active phase has product measure whereas nontrivial correlations are found numerically in the 
absorbing phase. It is argued that the system relaxes to the latter phase slower than exponentially. 
The critical behavior of this model is found to be different from that of the other known universality 
Ci ' classes. 
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I. INTRODUCTION 



Active-absorbing state transitions form an important class of nonequilibrium phase transitions whose best studied 
example is directed percolation (DP) Q. A critical behavior different from DP is possible in the presence of additional 
symmetries as in parity-conserving DP 2] . Sandpile models in which the activity occurs when the local density exceeds 
a certain threshold can also exhibit such transitions, provided the total density ("energy") is conserved Depending 
on the nature of the coupling between the activity and the local density, the conservation law may Q or may not [3 
change the critical behavior from DP. Further, the non-DP behavior can belong to different universality classes such 
(_J , as C-DP which describes stochastic sandpiles with linear coupling jfj, tj| . 

The models mentioned above have resisted an exact solution so far and have been studied extensively using numerical 
simulations and field-theoretic techniques. In this article, we introduce and study the steady state and the dynamics 
of a simple FESM for which some exact results can be found. As we shall see, the critical behavior of this model is 
different from that of the other known classes of FESMs. 

In Section m we define the model and study its steady state. In our model, only one particle hops out of a site 
£j ' to a nearest neighbor if the number of particles at that site exceed a certain threshold. In the high-density active 
phase, this model defines a zero range process (ZRP) whose steady state is known to be of product measure form 
in all dimensions H, 0. In the low-density absorbing phase, the steady state is not unique and depends on the 
initial conditions. We determine the critical behavior of this model by studying an order parameter and a correlation 
1 function close to the critical density. The time-dependent properties of this model are described in Section IIIII We 
' argue that the relaxation dynamics are the same as that of a two-species annihilation process; in particular, we find a 
stretched exponential decay in the inactive phase, which has been seen numerically in other sandpiles also |Toj . Our 
*0 ■ arguments are supported by Monte Carlo simulations. 
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II. THE MODEL AND ITS STEADY STATE 



The model is defined on a d-dimensional hypercubic lattice of length L with periodic boundary conditions. Each 
1 ■ site can hold an arbitrary number of particles with unit mass. We call a site active if the mass at this site is greater 
than mthr otherwise inactive. A particle hops out of an active site to a nearest neighbor with equal probability. In 
continuous time, the hopping rate u(m) at a site having m particles is given by 

u{m) = i l J 2d > m > m ^ (1) 
I , m < mthr ■ 
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The dynamics conserve the total density p — M/V where M is the total number of particles in the system and V = L d 
is the volume of the system. In one dimension with mthr — 1; this model can be mapped to a conserved lattice gas 
(CLG) model by regarding the sites as holes and mass as particle clusters. In the CLG language, a particle with one 
occupied neighbor can only hop to an empty nearest neighbor whereas the isolated particles do not move 

To see the phase transition, we define all the particles at a site to be immobile if the mass m < m t hr at this site. 
For m > mthr, the first mthr particles are said to be immobile and the rest m — m t hr mobile. If the initial number 
of particles is less than m t h r V , a mobile particle diffuses around till it reaches a site with m < mthr whereupon it 
becomes immobile. Thus, the number of mobile particles eventually becomes zero. On the other hand, if the total 
number of particles in the initial state exceed mthrV, then the number of mobile particles decrease till it reaches 
M — m t h r V. Thus, there is a phase transition at the critical density p c = m t hr from an inactive phase with only 
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immobile particles to an active one with a finite number of mobile particles, as the total density is increased. For sake 
of simplicity, we choose m t hr = 1 in the following discussion. 



A. The active phase 

We first note that the probability P/c(0, t) of having m = at site k at time t is a monotonically decreasing function 
in time since 



dPk(0,t) 



^il^w*' '" 1 ''' ' (2) 



dt 2d 

S m=2 

where Pk^k+sim' ' , m, t) is the joint probability that the site k and its neighbor k + S have mass ml and m respectively at 
time t. In the active phase, since p > 1, it follows that the probability P(0) of having an empty site is zero. It is then 
easy to check that the condition of detailed balance holds with the steady state distribution P(C = {mi, ...,my}) 
being equally likely Since the number of ways in which M — V particles can be partitioned in V cells is given by 
Z = (y-i) 0' we nave P(C) = 1/Z. Using a similar reasoning, the mass distribution P(m) = ( M y m ~ ls ) /Z can be 
obtained. Thus, the activity S = J2 m =2 p ( m ) = i M ~ V ) l( M ~ !)• Near tne critical point, S ~ (p- p c ) p with (3 = 1 
in all dimensions. Intuitively, as p — > p + 1/V, due to the absence of empty sites, the newly added particle typically 
becomes a mobile particle, leading to a linear growth of activity. This argument fails for the models studied in in 
which more than one particle can leave an active site thus creating empty sites. In such cases, one may expect the 
activity to grow nonlinearly. 

Since P(0) = for our model, the exponent f3 is expected to be unity even if u(m) is mass-dependent for m > mthr- 
This can be shown by realising that the active phase is a special case of ZRP for which P(C) ~ nLi /( m <0 where 
the marginal /(wfe) is given by 



/w = (n:;= 2 iM^) , m k >2 (3) 

yi , m k = U, 1 . 

The distribution P(m) = v m f(m)/Z is obtained by using the measure P(C) above with fixed density constraint. Here 
the normalisation constant Z = J2 m _i v m f(m) and the fugacity v is determined by the conserved particle number 
condition which can be written as 

p-l = |^(m-l)v m /M . (4) 

m=2 

Close to the critical density, the fugacity v — > and we can approximate p — 1 by the first term in the sum above. 
Using this in the expression for the activity S = J2 m =2 v m f(m)/Z for v — > 0, it follows that (3=1. 

Typically, the density-density correlation function C(r, p) = {mo m r ) — p 2 decays exponentially with r and can 
be used to define a static correlation length £ ~ (p — p c )~ l/± . Since this phase has product measure, the correlation 
function C(r,p) is zero in the thermodynamic limit for all r. This implies that the correlation length £> = and the 
exponent v±_ is undefined in this phase. 



B. The inactive phase and the critical point 

For p < p c , the system is in the inactive phase and the number of absorbing states is given by (£). Unlike the 
active phase, here the steady state is neither unique nor has product measure even if the initial state has. To see this, 
let us prepare an initial state in which a site is occupied by a monomer with probability pi , a dimer with probability 
P2 and remains empty with po — 1 — p\ — pi- Then, for p\ — 0, the final state will have one-clusters with only even 
number of ones whereas for p\ ^ 0, the configurations with odd number of ones are also allowed. Thus different 
initial conditions give different solutions for P(C) so that the steady state is not unique. Further, this is clearly not 
a product measure state for which the weight of all configurations at a fixed density is same. 

Using Monte Carlo simulations in one dimension, we measured the correlation function C(r,p) = (mo m r ) — p 2 
where m is either or 1. Since the steady state depends on the initial condition, the angular brackets denote spatial 
averaging for a given initial condition. We initially distribute particles independently at each site with the mass chosen 
from a Poisson distribution. As shown in Fig.^ the behavior of the correlation function C(r, p) when averaged over 
such initial conditions is similar to when only spatial averaging is done for a fixed initial configuration. Therefore, in 
the following discussion, we will carry out ensemble averaging as well. 
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The data in the inset of Fig.[2]for C(r, p) at various densities for a large, one-dimensional system supports the claim 
that the product measure does not hold in this phase. For densities away from unity, the correlation function decays 
fast and as one approaches the critical point, the decay gets slower. The magnitude of C(r, p) at small r drops for 
p — * 1 since at the critical density, we must obtain C(r,p c ) = 0. Also, since the total mass is conserved, we have 

(Mitiq) — Lp 2 — (toq) + (nio m i) ~\~ ••• + ( m o m L—i) ■ (5) 

Thus, the correlation function obeys the sum rule ^ r C(r, p) = due to which C(r, p) is not positive for all r. Finally, 
as shown in Fig. [21 a data collapse for C(r,p) for densities close to the critical density can be obtained if we assume 

CWJw^Cr/W > U = (l-P)- 2 ■ (6) 

The above scaling form is consistent with vanishing C{r,p) at the critical point. Thus, the static correlation length 
£< ~ (Pc ~ p)~ Ux diverges at the critical point with exponent j/f w 2. 



III. TIME-DEPENDENT BEHAVIOR 



Starting from an initial state in which particles are Poisson distributed, we study the temporal decay of the activity 
at and below the critical point. The activity S(t) at time t obeys the following equation 



dS{t) 



dt 2d 
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9d ^ 



,k+s (2, to, t) - Pk,k+s (2, 0, t) 



,m=2 



(7) 



The first term on the right hand side (RHS) represents the gain in the activity when a particle hops out of an active 
site to an inactive site, whereas the last two loss terms correspond to a particle leaving an active site with two particles. 
At low densities, since the probability of having large mass is exponentially small, the first two terms on the RHS of 
Eq.0 can be ignored and S(t) can be approximated by P(2, i). Denoting to = by A, m = 1 by (f> and m = 2 by B, 
the last term describes a two-species annihilation reaction A + B — ► <f> where A is a static and B is a diffusing species. 
A similar mapping to reaction-diffusion system for the models in Q involves more complicated reactions due to the 
creation of empty sites. In the following, we will denote the number of species A and B at time t by nA(t) and ns{t) 
respectively. 

At the critical density p c = 1, we expect that ns{0) ~ tia(0) and in a finite domain of length £, they differ by 
±0(i d l 2 ). Since the species B performs diffusive random walk, at large times in an infinite system, we have [T^ | 

w,*)~lrf 4 > d <t (8) 



J- 1 ,d>4 

where refers to the critical density in the thermodynamic limit. For a finite system of size L, the activity S at 
density p^ is expected to be of the scaling form, 

S(p?,t,L)*t- e H{t/L*) , (9) 

where the scaling function H{x) is a constant for i« 1 and grows as x e for x 3> 1 3- In other words, the activity 
typically decays as a power law in time to a system size dependent constant L~~ ze . As shown in Fig.Q, in one 
dimension, the activity S is of the scaling form in Eq.0 with z w 2 and the exponent 9 given by Eq.JSJ. However, 
the scaling function H(x) decays exponentially for x ^> 1 as seen by the linear decay of scaled activity on the semilog 
scale in Fig.®. In two dimensions, our simulations support the scaling behavior Eq.© with the scaling function 
decaying as a stretched exponential. The finite size scaling of the activity S oc L~ z0 in the steady state does not 
hold since the exact expression for S = (M — V)/(M — 1) gives p c = 1 for any L (unlike other sandpiles where p c is 
expected to be size-dependent). 

For p <C Pc, we expect rts(0) <C ^a(O) so that ua does not change appreciably due to the annihilation reaction and 
we are led to consider the problem of diffusing B species in the presence of static traps A. In this case, at large times, 
the density S(t) decays as a stretched exponential Q, 

S(t) ~exp[-(iA d ) Q ] , forp<p c , (10) 

where a = d/(d + 2) and i<j ~ (1 — p)~ 2 / d . This slow decay arises due to large trap-free regions which are rare 
but enhance survival significantly. The inset of Fig.® shows the decay of S(t) in one dimension in accordance 
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with Ea. i|10fl . This behavior of slow relaxation to the steady state is similar to that observed numerically in various 
stochastic FESMs in [To| . 

We next consider the effect of a perturbation in the steady state by considering the correlation function G(r, t) = 
(?7(0, 0) 77 (r, t)) where 77(7% t) is one if the site r is active at time t and zero otherwise. We perturb the steady state by 
constructing an initial condition with a single mobile particle placed at the origin. At the critical point, this particle 
executes a random walk and the correlation function G(r, t) is equal to the probability that the walker is at (r, t) so 
that 

G(p C) r,t)«^i^ex P (-r 2 /2t) . (11) 

In the absorbing phase, there is a finite density 1 — p of the vacant sites and the typical distance 7*0 between them 
scales as (1 — p)~ x l d . Therefore, the single mobile particle at the origin diffuses for t <C r§ so that the correlation 
function G(f,t) ~ e~ r ~/ 2t ,r < r increases with time. However, for t ^> r%, the mobile particle gets trapped and 
G(f,t) decays as in Ea. lfUT|l . Usually, the auto-correlation function G(0,i) is expected to decay exponentially, i.e. 
G(0,i) ~ e~*/ T where r ~ (p c — p)~ VH and the exponent v\\ obeys the scaling relation 2/11 = v±z 0. The stretched 
exponential behavior of G(0, t) in our model implies that is infinite. Alternatively, one can define this exponent 

via td ~ {pc — p) with vf? = 2/d. Since, in our model, z « 2 and the correlation length exponent 2/^ » 2 in one 
dimension, either definition of 2/5 is inconsistent with the above scaling relation and we conclude that it does not hold 
in the absorbing phase. 



IV. CONCLUSION 



We introduced a sandpile model whose simplicity allowed us to determine various static and dynamic exponents 
exactly. These exponents differ from those of the C-DP class for which, in one dimension, f3 ~ 0.29, 2/^ rj 1.33, 
z « 1.55 and 9 m 0.14 0. We considered the correlation lengths and defined above and below the critical 
density respectively and found that while is zero, diverges as p — > p c — . However, a further detailed study of 
the nontrivial correlations present in the inactive phase would be interesting. We also studied the temporal properties 
of this system by relating it to a well studied reaction-diffusion system via an argument similar to that used in [l5| 
to study avalanche size distribution in a sandpile model. 

One can also consider the biased version of the above model in which a particle moves preferentially to a nearest 
neighbor. The steady state of this model is the same as that for the symmetric case discussed above. For the 
asymmetric case, using the known results for the two-species annihilation problem with B species drifting with a 
nonzero speed, we expect that at the critical point, S(t) decays as a power law with an exponent d/2 for d < 2 and 
exponentially in the subcritical regime |16|. Our numerical results are consistent with the argument. 
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FIG. 1: Plot of C(r,p) vs. r for p < p c in one dimension with Poisson distributed initial mass. The data shown with points 
corresponds to different initial seeds for the random number generator. The data averaged over 1000 initial conditions is shown 
with broken line. 
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FIG. 2: Inset: Plot of C(r,p) vs. r in one dimension for p < p c with L = 16384 to show that the product measure does not 
hold in the inactive phase. Main: Data collapse for the scaled correlation function £ C(r,p) vs. r/£ with £ = (1 — p)~ 2 . 
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FIG. 3: Data collapse for the scaled activity t 1 ^ 4 S(t, L) vs. t/L 2 at p — p c in Id in accordance with Eq.Q. The inset shows 
the temporal decay of S(t) for p < p c as in Eg. llOt . 



